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ABSTRACT 

We perform a local (short-wavelength) linear stability analysis of an axisymmetric col- 
umn of magnetized plasma with a nearly toroidal magnetic field and a smooth poloidal 
velocity shear by perturbing the equations of relativistic magnetohydrodynamics. We 
identify two types of unstable modes, which we call 'exponential' and 'overstable', re- 
spectively. The exponential modes are present in the static equilibria and their growth 
rates decrease with increasing velocity shear. The overstable modes are driven by the 
effects of velocity shear and dominate the exponential modes for sufficiently high shear 
values. We argue that these local instabilities can provide an important energy dissipa- 
tion mechanism in astrophysical relativistic jets. Strong co-moving velocity shear arises 
naturally in the magnetic acceleration mechanism, therefore it may play a crucial role 
in converting Poynting-flux-dominatcd jets into matter-dominated jets, regulating the 
global acceleration and collimation processes, and producing the observed emission of 
blazars and gamma-ray bursts. 
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1 INTRODUCTION 

Many high-energy astrophysical phenomena are powered by 
jets that form in the vicinity of compact objects and are 
accelerated to relativistic speeds (blazars and gamma-ray 
bursts (GRBs)) or very high subrelativistic speeds (X-ray bi- 
naries, pulsar wind nebulae). The currently standard model 
of magnetic acceleration requires that jets be initially dom- 
inated strongly by their Poynting flux, and that accelera- 
tion takes a few orders of magnitude in distance scale (e.g., 
Komissarov et al. 2007; Lyubarsky 2009). In this process, 
the jet expands laterally far beyond the light cylinder of the 
central object, stretching the magnetic field into a nearly 
toroidal configuration. Magnetically-dominated plasma with 
a toroidal magnetic field is known to be highly suscep- 
tible to Current-Driven Instability (CDI). These instabili- 
ties cause severe problems in laboratory experiments on 'Z- 
pinch' plasma systems. However, in the context of astrophys- 
ical jets they provide both a challenge and an opportunity. 
The challenge is not to destroy the jet as it propagates over 
several orders of magnitude in distance. The opportunity is 
to provide a mechanism of energy dissipation — by tapping 
free magnetic energy and converting it to particle energy — 
that can power high-energy emission from blazars and GRBs 
(Giannios & Spruit 2006). 
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Global, three-dimensional general-relativistic magneto- 
hydrodynamical (3D GRMHD) numerical simulations of 
jets, tracing them from their formation around spinning 
black holes over distances of thousands of gravitational radii, 
indicate that they are marginally stable to long-wavelength 
modes of CDI, including the most dangerous pinch (m = 0) 
and kink (m = 1) modes (McKinney & Blandford 2009). On 
the other hand, it is expected that the growth rate of CDI is 
similar to the Affven crossing time scale of the jet, and thus 
should not be sensitive to the wavelength. Short- wavelength 
(local) modes are particularly important for the problem of 
energy dissipation, since the dissipation rate scales like the 
inverse square of the wavelength. These modes were pre- 
dicted analytically to play an important role in AGN jets 
and pulsar wind nebulae (Begelman 1998), and recent local 
numerical simulations have confirmed their general proper- 
ties (Mizuno et al. 2009, 2011a; O'Neill et al. 2012). 

Most analytical studies of CDI have been performed 
in the force-free approximation, i.e. neglecting the pressure 
and inertia of the matter content of jet plasma (Istomin 
& Pariev 1996; Lyubarskii 1999; Tomimatsu et al. 2001; 
Narayan et al. 2009). This approximation is valid in the in- 
nermost regions of jets, within the fast magnetosonic point, 
at which the plasma magnetization a is still substantial. Be- 
yond the fast magnetosonic point, the magnetic nozzle effect 
results in further gradual jet acceleration, and a decrease in 
the jet magnetization to the order of unity (Begelman & Li 
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1994). In this jet region, full MHD force balance is necessary 
to study plasma stability, as was done by Begelman (1998). 
However, that study assumed a static background configura- 
tion, while it is known that the efficiency of jet acceleration 
is strong function of cylindrical radius (Komissarov et al. 
2007; Tchckhovskoy et al. 2008). Within the jet accelera- 
tion region, a strong shear of the poloidal velocity can be 
expected to arise, and this effect needs to be incorporated 
in the local stability analysis. 

Velocity shear is the driving factor of the Kelvin- 
Helmholtz Instability (KHI) , which has been studied exten- 
sively in the context of astrophysical jets. Most analytical 
studies of KHI have adopted a vortex-sheet background con- 
figuration, i.e. a parallel discontinuity of the velocity field 
(Turland & Scheuer 1976; Blandford & Pringle 1976; Ferrari 
et al. 1978; Birkinshaw 1984; Hardee 2007). In some works, a 
more gradual transition region of finite width was considered 
(Ray 1982; Ferrari et al. 1982; Birkinshaw 1991). However, 
we have found no analytical study that deals with a strictly 
smooth velocity shear, which is natural for the local stability 
problem that we study in the context of CDI. Several nu- 
merical studies have investigated the effect of velocity shear 
on CDI or the interaction between CDI and KHI modes. 
Some of them have adopted the vortex-sheet approximation 
(Mizuno et al. 2011b), while some use a gradual velocity 
shear with a well-localized region of strong shear generating 
typical KHI modes (Baty & Keppens 2002). For moderate 
shear values, the distinction between CDI and KHI modes 
becomes quite vague. 

In this work we study local CDI modes of cylindrical 
MHD equilibria in the presence of strong but gradual ve- 
locity shear. In Section 2, we derive the dispersion relation 
for a sheared plasma column, which is given by Equation 
(25). In Section 3, we explore the complex space of mode 
frequencies, looking for the fastest-growing mode as a func- 
tion of the background equilibrium parameters. In Section 
4, we discuss the relevance of our results to global models of 
astrophysical jets. We summarize our results in Section 5. 

In this work, we use cgs units with c = 1. In the linear 
analysis, all quantities describing the background configura- 
tion are denoted with subscript '0' (zero-order terms), while 
all quantities describing the linear perturbation are denoted 
with subscript '1' (linear terms). Vector quantities are de- 
noted with a bold font. 



Bo,^ is radially dependent, with logarithmic derivative — 
<9inr(mi?o,0)- The poloidal magnetic field Bo,z <C B ^ is 
a constant 1 . There is no background radial magnetic field 
(5o,r = 0). The gas pressure po is in equilibrium with the 
magnetic field: 



Cpo = -(a<t> + 1 



47T 



4tt 



(1) 



where £ = <9i nr (lnpo)- Here, we keep the term including vo,z 
as it will be eventually differentiated with respect to r. 

We introduce linear velocity perturbations of the form 



v\ (t, r, (f>, z) — vi exp(wirf + ilr + imtf) + ikz) 



(2) 



These induce perturbations of the gas pressure (pi) and the 
magnetic field (Bi). The pressure perturbation is found by 
perturbing the equation of state dt(pp -7 ) = in combina- 
tion with the continuity equation dtp + V • (pv) = 0: 

d t pi + {v! ■ V)po + (v Q ■ V)pi + 7po(V -vi) = 0, (3) 

the solution to which is 

= t(C + 7K r -m 7 ^-7«i,irkz (4) 
(w + kvo,2.)r 

where ui,i r kz = (lvi, T + kvi, z )r is the perturbed velocity com- 
ponent that we will use instead of vi, x . The magnetic field 
perturbation is found by perturbing the Maxwell equation 
V x E — —dtB together with E = B x v. 

d t Bi = (Bi ■ V)v + (-Bo ■ V)wi - (wo ■ V)Bi 

-(wi-V)B -(V-wi)B . (5) 
The solution is: 



(m + T?)ttl, r 
Dl,r = 7 ~T ~ t >0,(j> 



Bi,4 



ia<t>Vi,i + (yvi,4> — m.irkz) 
{u> + kvo, z )r 



B, 



B _ i[(m + ri)ilr + rj\vi, T - m(j)V\^ - «i,i rk z) B 
" (u> + kv , z )kr 2 "' 4 

i(m + r/)uvi^ 



(ui + fcuo,z) 2 r 



(6) 
(7) 

(8) 



where rj = kr(Bo,z/B ^). We note that _Bi jZ contains a 
non-vanishing term proportional to the background veloc- 
ity shear. 



2 DERIVATION OF THE DISPERSION 
RELATION 

Here, we perform a rigorous local linear perturbation analy- 
sis of a cylindrically symmetric plasma column filled with a 
predominantly toroidal magnetic field, in equilibrium with 
a gas pressure gradient in the radial direction. We gener- 
alize the results of Begelman (1998), taking into account a 
smooth poloidal velocity shear. 

We introduce cylindrical coordinates (r, <j>, z). Calcula- 
tions are made at a certain radius r in the local co-moving 
frame. Hence, the local poloidal plasma velocity is vo,z = 0, 
while the local poloidal velocity shear is u — £h nr vo,z > 
(as we show later, the choice of the sign of u is not sig- 
nificant). We assume no background rotation (vq^ = 0) 
or radial expansion (vo,i = 0). The toroidal magnetic field 



2.1 Energy- momentum conservation law 

Energy-momentum conservation is governed by the equation 
d„T^ = 0, where 

T tm = T 2 w-p+ — <E 2 + B 2 ) (9) 
T 01 = T 2 wv x + ^-(Ex BY (10) 

47T 

T ij = r 2 uwV + (p + e2 ^ b2 ^ s 13 

-±-.(e 1 E 3 +B l B 3 ^j (11) 

1 Introducing a radial dependence of the poloidal magnetic field 
has no effect on the dispersion relation under the ordering of terms 
described later on. 
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is the stress tensor in the ideal relativistic MHD regime, 
T = (1 — v 2 )^ 1 ^ 2 is the Lorentz factor and w is the relativistic 
enthalpy. The combined energy-momentum equation is 



Pro (d t + v ■ V) v + (vd t + V)p 

-p e E-jxB + (E-j)v = 0, 



(12) 



where p c is the electric charge density and j is the elec- 
tric current density, which can be found from the Maxwell 
equations: 



±(VxB-d t E). 



(13) 
(14) 



The perturbed energy-momentum equation is 
ToWo [dtvi + (vi ■ V)«o + (vo ■ V)wi] 

+V Q d t pi + Vpi - p e ,lE - Pe.O-Bi 

-ji x B„ - j x B-l + {Ei ■ jo)v + (E ■ j^vo = . (15) 

We restrict this analysis to perturbations of short wave- 
length in the z direction, i.e. kr S> 1. At the same time, we 
assume that Bo,z -C Bo,f so that ij ~ 0(1). For the wave- 
length in the r direction, we assume only that I < O(k). 
We further assume that m,u ~ 0(1); and introduce the di- 
mensionless frequency uj = ojr ~ 0(1), the magnetization 
parameter oo = Bq/(Attwo) ~ 0(1), and the sound speed 
u s = JPo/wo ~ C(l)- Wc can now calculate each term of 
the energy-momentum equation (see Appendix A). 

The azimuthal momentum equation leads to (keeping 
the leading terms only): 

(mv 2 + 770-0) t>i,r + (m 2 v 2 + r/ 2 a - Co 2 ) ivi^ 

+ (mv 2 - 770-0) ivi.irkz = , (16) 
while the poloidal momentum equation gives: 
(« s 2 - 0- o )vi,r + (mv 2 - r/a )iv lt<t> 

+ (v 2 + <To)iVl,lrkz = 0. (17) 

We can combine these two equations to derive explicit for- 
mulas relating the perturbation velocity components: 

2(m + ri)v 2 a 



Vl } <h = 



«l,lrkz = 



(m + rj) 2 v 2 a Q - Co 2 (v 2 + <j 
(;rj 2 - m 2 )v 2 a - (v 2 - a»)C 



■Wi,r , 



-IVl,, 



(18) 



(19) 



(rn + ri) 2 v 2 a - Cb 2 (v 2 + <t ) 

This means that the components «i jr , Vi^, Vi.irkz of the per- 
turbation velocity are of the same order. We can now com- 
pare the orders of terms involving different velocity compo- 
nents. 

To 0(vi), the radial momentum equation multiplied by 
a factor (r/wo) is 



din rPl 
w 



+ 



/1 \~ il 
(1 + c )uj — -a « 



ivi, T — i(m + 17)0-0 



£1, 
B, 



0,4, 



(20) 



+K + 2)ao _ +a ^__ + ___j =0 . 

In order to eliminate the pressure term, we calculate the 
r derivative of the poloidal momentum equation (Equation 
15) multiplied by a factor (ir/kwo) to 0(vi): 



din rPl 
Wo 



+ 



I 2 ,, il 



IVl, 



-a.4,00- 



Bi, 



o 



ftnr-61, 



m dlnrBl, 



0. 



Bo,<j> V Bo,4> kr Bq^ 

Summing Equation (21) with Equation (20), we obtain 



(21) 



I 2 \ 

1 + ^2 I (1 + ffo)iw?Ji,r - i(m + 77)0-0 



+2o 



B14 



(m + 77) 
H ; 00 



-Bo,0 

dlnrB\^ 



(22) 



0. 



We need to calculate <9i nr i3i,z, but only to O(krviB ,4,): 
di nr Bi >z 



B, 



o,4> 



— kr(m + 77) 



I 2 

'Y 2+2 



il Co\ u 
k u j Co 1 



^.(23) 

Co 



After substituting the perturbed magnetic field to Equation 
(22), we obtain: 



+2o-„ 



1 + ) + - ( m + "ifvo] 



+ (rn + 77) 1 



Wi,, 
Co 



(24) 



+2o 



il Co^ 
k u J LO z 



Cb 

"l.lrkz) 



0. 



In the last term, we substitute the perturbation velocities 
using Equations (18) and (19). Finally, dividing Equation 
(24) by a factor (2icroVi ;I /oj), we obtain the dispersion rela- 
tion: 



K 1+ 



1 + — ) Cb z — (m + tjY 



(To 



+ (m + 77)' 
(m + r)) 2 v 2 a 



^ il Co 
k u 
Co 2 '(ao - 



(m + r)) 2 v 2 a - Co 2 (a + v 2 ) 







(25) 



(see Appendix B for a notation look-up table). For u = 0, 
the above equation corresponds exactly to Equation (3.32) 
of Begelman (1998). The velocity shear term can be traced 
back to the extra perturbed poloidal magnetic field _Bi jZ (see 
Equations 8 and 23) via the j\ x Bo force (Equation A6). 
The (rn + rj) 2 factor indicates that it contributes only for 
rn 7^ in the presence of a purely toroidal field (77 = 0), or 
requires a finite background poloidal magnetic field (77 ^ 0). 
The velocity shear terms become dominant for lo — > 0. 

The KHI modes, driven by gas pressure in the presence 
of velocity shear, are suppressed in the short wavelength 
approximation. This can be seen by noting that the vod t pi 
force would contribute a term —coupi/kwo to Equation (21), 
which is of order 0(v\/kr). Even in the absence of the mag- 
netic fields, it would be dominated by the acceleration term 
by one order in kr. The KHI modes would be important only 
when |u| > 0(kr). 



3 STABILITY ANALYSIS 

In the general case, the dispersion relation derived in the 
previous section (Equation 25) is equivalent to a sixth-order 
polynomial over the complex space of Co. Hence, we cannot 
perform a fully analytic stability analysis, as is possible for 
the case of a static column, when Equation (25) reduces to 
a quadratic equation for Co 2 with only real solutions (Begel- 
man 1998). Substituting Co — Cor — iCoj, where wr,wi G 1Z, 
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we can separate Equation (25) into its real and imaginary 
components (we also substitute a = (m + r)) 2 and d — l/k, 
see Appendix B): 



(1 + rf 2 



! + -)(£& 



~2x 

• Wi ) — a 



(ui R — aij ) aduu)i 



(to 



,7^2 



+ 



/ - 2 



- 2 
- Wi 



2 /~2 . ~2\2 4 | ,1-2-2 2 



[aw|cr - (a) 2 , - cDj 2 )(cro + w|)] 2 + 4w R w 2 (cr + «f ) 2 



(1 + d 2 



1 + ^ 

O"0 



_ . 2au iDrWi 

WrWi — „ _„ 

(<4 + w i ) 

4a«gCToWR(Di 



\avioo 



■w 2 )(ao + « 2 )] 2 +4w 2 a; 2 (ao+« 2 ) 2 



(26) 
= 0, 

(27) 
= 0, 



respectively. The imaginary equation has a trivial solution 
at lob. = 0. Hence, for purely exponential solutions it is suf- 
ficient to solve only Equation (26). We further note that 
Equation (26), and Equation (27) divided by cD R , contain 
only quadratic terms in £ir. Hence, it is sufficient to look for 
solutions with lob > 0. Also, we are interested solely in so- 
lutions with uj\ > 0. In the case of vanishing velocity shear, 
the left-hand side of Equation (27) is always positive, hence 
there are no solutions with wr > 0. 

In Figure 1, we show typical solutions to the system of 
Equations (26) and (27). We fix the values for a^, u, ao and 
v s , and for selected values of a we continuously vary param- 
eter d. For all values of a and d, there are some overstable 
modes (with Cj-r, 0), which form two separate branches in 
the (wr,(I>i) plane. However, for sufficiently small values of 
a, there exists a class of exponential solutions, which dom- 
inate the overstable solutions in terms of the growth rate 

til. 

In Figure 2, we show the fastest-growing modes as a 
function of a, selected by varying d, for a fixed value of 
a,£ = —0.5 and several values of the velocity shear u. For 
u = 0, such a system is known to be unstable only for 
modes with a < 1, which is consistent with the stability cri- 
terion a < 2(a^> + 1) (see Equation (4.2) in Begelman 1998). 
These unstable modes are purely exponential. For small non- 
zero values of u, overstable modes arise for a > oi, where 
di < 1. These modes are slower than the exponential ones 
for u = 0.1, but of comparable growth rate for u — 0.2. For 
u > 0.2, the overstable modes dominate the exponential so- 
lutions, and their growth rate is roughly proportional to u. 
The value of a corresponding to the fastest-growing modes 
increases with increasing u for the overstable solutions, and 
decreases with increasing u for the exponential solutions. 
The overstable solutions generally have their oscillatory fre- 
quency cDr > wi. 

In Figure 3 we show the fastest-growing modes selected 
by varying both d and a, as a function of the velocity shear 
u. In the case of = —0.5, it is confirmed that the fastest- 
growing modes are purely exponential for u < 0.2 and over- 
stable for u > 0.2. We confirm the behavior of a with u 
identified in Figure 2. In the case of = 0, the exponential 
modes with a = d = strongly dominate over the overstable 
modes and the velocity shear has very little impact on the 
former. The case of — —1 corresponds to marginal sta- 
bility in the absence of velocity shear. In this case, there are 
no purely exponential solutions, but the system is prone to 




Figure 1. Sample solutions to the system of dispersion relations 
given by Equations (26) and (27). For each value of a = (m + »?) 2 , 
we plot solutions for a wide range of d = l/k. Both panels show 
the same set of solutions in different projections. 



overstable modes. Both the growth rate and the oscillation 
rate are proportional to the velocity shear, with approxi- 
mate scaling relations u>i ~ 0.62-u and ur ~ u. Overstable 
modes are present also for — — 2. 

Figure 4 shows the dependence of the fastest growth 
rates on the magnetization ao and the sound speed v s . When 
changing the plasma magnetization, we need to adjust the 
toroidal magnetic field distribution to maintain a realistic 
radial force balance. We use Equation (1), setting vo, z = 
and po = wo/4 to obtain a relation (aj, + 1) = — £/(4<ro)- 
For £ = 0, i.e. a uniform distribution of the gas pressure, we 
have = —1. In Figure 4, we show a non-trivial case of 
C = — 2. We find that the plasma magnetization has a signifi- 
cant effect on the growth rates of the exponential modes and 
the minimum velocity shear for which the overstable modes 
dominate the exponential modes, but it has a weak effect 
on the growth rates of the overstable modes. For a given 
value of the velocity shear u, the growth rate decreases with 
increasing ao, i.e., highly magnetized plasmas are more sta- 
ble. The limit of ao ^> 1 is called the force-free limit. In 
this case, the dimensionless growth rate for the exponential 
modes tends to zero, and the overstable modes can domi- 
nate the exponential modes for very small values of u. For 
moderate values of ao, the growth rate is sensitive to the 
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100 



a = (m + n) 



Figure 2. Fastest-growing modes as functions of the parameter 
a, and for different values of the velocity shear u, obtained by 
varying the parameter d. Thick lines show the growth rate &i, 
while thin lines show the oscillation frequency o)r. The dashed 
vertical lines mark the transition between the exponential and 
overstablc modes, with overstable modes dominating to the right 
of the line. 



sound speed (i.e. to the plasma temperature), being higher 
for colder plasmas. 

Our results indicate that the overstable modes can be 
driven for any value of a^, ao and v s , and that they domi- 
nate other solutions for a sufficiently large value of u. Both 
the growth rate of the overstable modes and the oscillation 
frequency are roughly proportional to the velocity shear. 



0.8 
0.7 
0.6 
0.5 
iS 0.4 
0.3 
0.2 
0.1 



c = 1 
v c 2 = 0.333 




Figure 3. Fastest-growing modes selected by varying parameters 
d and a, as functions of the velocity shear u for different values 
of parameter ct^. Top panel: the growth rate Co\. Middle panel: 
the oscillation frequency Bottom panel: the optimal value 
of parameter a. For the cases having exponential solutions, the 
dashed vertical lines mark the transition between the exponential 
and overstable modes. 



4 DISCUSSION 

In order to evaluate the importance of the velocity shear in 
affecting the stability of astrophysical jets, we first need to 
estimate the magnitude of this parameter. Our calculations 
in previous sections were performed in the local co-moving 
frame (in this section, we denote the quantities calculated in 
the co-moving frame with a prime). The co- moving velocity 
shear is directly related to the transverse gradient of the 
Lorentz factor calculated in the external frame: 

' / ^2 a 1 ainr 

u = ra r v . = T rd r «o,z = tt, • (28) 

«o,z dlnr 

In relativistic jets, a co-moving velocity shear value of u' = 1 
corresponds to the doubling of the Lorentz factor (as mea- 
sured by an outside observer) with a doubling of the radius. 
The co-moving velocity shear is significant when the Lorentz 
factor varies by a sizable factor across the jet. 

Relativistic jets are thought to be accelerated by the 
action of magnetic stresses in an extended region dominated 
by the Poynting flux. Both analytical studies (Vlahakis & 
Konigl 2004; Lyubarsky 2009) and numerical simulations 
(Komissarov et al. 2007; Tchekhovskoy et al. 2008; Porth & 
Fendt 2010) show that this process does not operate with 
uniform efficiency across the jet radius. This results in a 





1 1 1 1 1 1 1 1 1 1 — 7-7? 1 

* J?' 













\<y'' o = 1/3, 0^ = 0.5 




°o = 1 > a o = "°- 5 - 
c = 10, Ofy = -0.95 




a = 100, o^ = -0.995 



0.2 0.4 0.6 0.8 1 1.2 1.4 
u 

Figure 4. Dependence of the fastest growth rates on the plasma 
magnetization 00 and the sound speed v s . The toroidal magnetic 
field distribution is chosen to satisfy {a^ + 1) = 1/(2<to)- The 
black lines show the results for v% = 1/3, and the gray lines show 
the results for = 0.1. 
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Figure 5. Top panel: three functions for the Lorentz factor T dis- 
tribution across a jet, motivated by numerical results of Komis- 
sarov et al. (2007) and Tchckhovskoy et al. (2008). Bottom panel: 
the corresponding distribution of the co-moving velocity shear u' . 



strong dependence of the Lorentz factor within the jet on 
the radial coordinate (equivalently on the polar angle). To 
illustrate these dependencies, in the top panel of Figure 5, 
we plot three models for the radial variation of the Lorentz 
factor, qualitatively based on the results of Komissarov et al. 
(2007) (models CI and C2) and Tchekhovskoy et al. (2008) 
(Appendix 4.1). In the bottom panel of Figure 5, we plot the 
corresponding values of the co-moving velocity shear given 
by Equation (28). We find that velocity shear values of order 
|w'| ~ 0.5 and higher can be expected in relativistic jets. This 
indicates that the effect of the poloidal velocity shear on the 
stability of jets is very important. A gradual jet acceleration 
also means that the importance of velocity shear can be 
expected to increase with distance. 

In our local approximation, the sign of the velocity shear 
has no effect on the maximum growth rate for a given sys- 
tem; only the magnitude of the shear is important. This is 
because in Equation (25), the linear term in u is coupled with 
the l/k ratio, and a change in the sign of u can be countered 
by a change in the sign of I (inward or outward propagation) 
or k (upstream or downstream propagation). On the other 
hand, Narayan et al. (2009) found that for global modes in 
the force-free approximation, the sign of the poloidal veloc- 
ity shear is important. They found that jets with positive 
velocity shear (Lorentz factor increasing with radius) are 
stable, while those in which the velocity shear changes sign 
at an intermediate radius, as in the example shown in our 
Figure 5 with a dotted line, are unstable. This fundamental 
difference in conclusions reached via global and local anal- 
ysis is not surprising. It demonstrates that even in globally 
stable jets, vigorous dissipation may result from the presence 
of local modes. 

The local stability criterion formulated by Begelman 
(1998) for a static plasma column, which can be written as 
a > 2(a^ + 1), is no longer valid if poloidal velocity shear 
is present. This is clearly illustrated in Figure 2. In the case 
of = —0.5, the stability criterion is a > 1. However, even 
for a very small value of u, there exist overstable modes with 
a > 1. Even for < —1, in which case the static equilibria 
are always stable, we find overstable solutions. 



A mode with a co-moving dimensionless growth rate of 
Cb'i will double its amplitude over the co-moving time scale 
of t' 2 — Q.7r/uj[. Over this time, an element located initially 
at z will travel a distance 22 in the external frame given 
by (22/2) ~ O.7(r0)/w{, where T » 1 is the local Lorentz 
factor and 6 = r/z is the polar angle. Since in relativistic 
jets we expect TO < 1, for Cj[ ~ 0(1) the growth length 
measured in the external frame is not particularly short, 
with (22/2) ~ O(l), with the exception of the core region 
close to the jet axis. For a typical co-moving growth rate of 
uj'i ~ 0.4, the energy dissipation associated with the local 
CDI can be expected to proceed over at least one order of 
magnitude in distance along the jet. A major caveat is that 
the growth of local CDI modes competes with lateral jet 
expansion. The question of whether jet expansion is able 
to quench the CDI modes must be addressed by a careful 
analysis of realistic jet models. 

The results of this work need to be confirmed by nu- 
merical simulations. We are currently preparing local sim- 
ulations of cylindrical plasma columns similar to those in 
O'Neill et al. (2012), adding a background poloidal veloc- 
ity shear. However, this issue also needs to be addressed 
by global simulations, like those in McKinney & Blandford 
(2009), but with a much higher resolution. 

Our study of the local development of CDI is still very 
simplified, as we do not take into account several effects that 
could potentially have an impact on the CDI growth rate. 
Some of these effects, like rotation or lateral expansion, are 
not supposed to be important at the scales we are inter- 
ested in. However, centrifugal forces due to jet collimation, 
enhanced by relativistic effects, can significantly change the 
radial force balance (Narayan et al. 2009) and should be 
included in future work. 



5 SUMMARY 

We have calculated the local dispersion relation (Equation 
25) for a cylindrical plasma column with smooth poloidal 
velocity shear (defined as u — rd r vo,z), extending the orig- 
inal results of Begelman (1998). The dispersion relation is 
a sixth-order complex equation for dimensionless frequency 
u = wr = wr — iui, which allows two types of solutions: 
'exponential' solutions, which have a purely imaginary fre- 
quency uii > 0; and 'overstable' solutions that in addi- 
tion have a non-zero real frequency wr 7^ 0. We explore 
the parameter space of mode wavenumbers and identify the 
fastest-growing modes for a given equilibrium. The growth 
rate of the exponential modes decreases with increasing ve- 
locity shear, while that of the overstable modes is roughly 
proportional to the shear. For a sufficiently high shear the 
overstable modes dominate the exponential modes. The 
transition between the two regimes depends on the distri- 
bution of the toroidal magnetic field (parametrized by — 
<9i nr (lni?o,0)), the plasma magnetization (<r = Bq/(4tywo)) 
and the plasma temperature (measured by the sound speed 
given by = jpo/wo)- In the force-free limit (ao 1), as 
a.$ ~ —1, the exponential modes become unimportant, and 
the overstable modes dominate over a wide range of shear 
values. 

These results can be applied to relativistic jets of AGN 
and GRBs, in the initial force-free region, in the main accel- 



© 0000 RAS, MNRAS 000, 000-000 



CDI with poloidal velocity shear 7 



eration region, and beyond. Studies of magnetic acceleration 
of jets predict a strong dependence of the jet Lorentz factor 
on the cylindrical radius, which translates to significant val- 
ues of the co-moving velocity shear. Hence, in many realis- 
tic jet configurations, the overstable modes identified in this 
work may dominate locally. Since these modes are present 
for every set of equilibrium parameters, they are very hard 
to suppress unless the underlying velocity shear is reduced to 
insignificant values. Local CDI may play an important role 
in regulating jet acceleration and collimation. In the mean- 
time, a significant fraction of the relative kinetic energy may 
be dissipated, transferred into ultra-relativistic particles and 
converted into non-thermal radiation that dominates the 
spectra of blazars and GRBs. 
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APPENDIX A: INDIVIDUAL TERMS OF THE 
PERTURBED ENERGY-MOMENTUM 
EQUATION 

Here, we calculate individual terms of Equation (15). We 
retain the terms linear in i>o, z only in the z component, as 
it will be differentiated with respect to r. The acceleration 
term is: 



r w [d t vi + (vi ■ V)v () + (v Q ■ V)vi] ~ 

wo , . , * „ 

~ ILdWoVl H (!fcrtJ[) >z Ul iZ + UVl^) z . 

Gas pressure forces are: 
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_ _ impi - 

Vpi = OrPir H <p + ikp\z . 

r 

Electromagnetic forces are: 
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APPENDIX B: NOTATION 



OJ — COB. — tUJl 


mode frequency 


CO = LOT 




\ . T7~l] k 


wflyp mini nprQ 
Wave 11LL111LJCIS 


A 111. 

a — 1 1 k 




Bo, 4, , Bo,z 


. 

magnetic field 


a<f, — Olnr(mi!o,</)J 




^ — U r ( R„ / R„ , \ 
77 — KT (J3Q,z/ £>0,4> ) 




a = (m + rj) 




Wo 


bpeClllC CIlLIldipy 


a = Bo/(4nw ) 


magnetization 


Po 


gas pressure 


C = 9lnr(lnp ) 




7 


adiabatic index 


vl = 7Po/w 


sound speed 




poloidal velocity 


W = <9l nr U0,z 


velocity shear 
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